Nonlinear saturation of electrostatic waves: 
mobile ions modify trapping scaling 



John David Crawford and Anandhan Jayaraman 
Department of Physics and Astronomy 
University of Pittsburgh 
Pittsburgh, Pennsylvania 15260 
(June 24, 1996) 

Abstract 

The amplitude equation for an unstable electrostatic wave in a multi-species 
Vlasov plasma has been derived. The dynamics of the mode amplitude p(t) 
is studied using an expansion in p; in particular, in the limit 7 — ► + , the 
singularities in the expansion coefficients are analyzed to predict the asymp- 
totic dependence of the electric field on the linear growth rate 7. Generically 
l-Efcl ~ 7 5 ^ 2 , as 7 —* + , but in the limit of infinite ion mass or for instabili- 
ties in reflection-symmetric systems due to real eigenvalues the more familiar 
trapping scaling \Ek\ ~ j 2 is predicted. 
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The evolution of an unstable electrostatic mode is a fundamental problem in collisionless 
plasma theory, and is perhaps the simplest nonlinear problem requiring a self-consistent 
treatment of the resonant interaction between waves and particles. When resonant particles 
interact with a large amplitude wave, then much of the behavior can be understood by 
analyzing the particle motion as if the wave amplitude were constant; this approximation 
linearizes the problem. [||] Alternatively, if the wave amplitude is sufficiently small then 
the initial instability can be predicted treated by neglecting the effect of the wave on the 
particles; this leads to conventional linear Vlasov theory. However, to describe the dynamics 
of the unstable mode which develops from a small initial amplitude into a final nonlinear 
state requires an analysis of the self-consistent and nonlinear interaction between the wave 
and the resonant particles. 

From a dynamical systems viewpoint even the simplest examples of instabilities in a 
Vlasov plasma have many unusual features related to the Hamiltonian character of the 
dynamics and the central role played by the neutrally stable continuous spectrum in the ap- 
pearance of the unstable modes. 0] These novelties are not present solely in Vlasov theory; 
entirely analogous features arise in models of unstable inviscid shear flows, in stability calcu- 
lations for certain classes of solitons, and in theories of large systems of coupled oscillators. 
||- || In the better understood setting of dissipative systems, the nonlinear evolution of 
the mode amplitudes can be described using an expansion in the amplitude of the unstable 
modes. Since the growth rates are very small near onset, nonlinear effects often act to satu- 
rate the instability before the amplitudes grow appreciably; for this reason such expansions 
have proved a powerful tool for studying the nonlinear states emerging from the bifurcation. 

It has long been hoped that similar methods could be applied to the Vlasov equation 
despite the absence of dissipation. However, for many years, efforts to construct such ex- 
pansions, even for the case of a single unstable mode, have been plagued by the fact that the 
nonlinear terms involved divergent integrals due to resonant denominators. Thus the calcu- 
lations appeared to breakdown precisely in the regime where the amplitudes of the unstable 
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waves were extremely small. In addition, efforts to regularize the the expansion coefficients 
inevitably led to theories that predicted scaling behavior for the saturated amplitudes that 
contradicted numerical results, fl2"|| - flEfl More precisely, these theories predicted that the 
electric field of the saturated mode would satisfy E ~ y^y as 7 — ► + , whereas numerical 
simulations find the exponential growth of the mode halted at an amplitude characterized 
by the "trapping scaling" E ~ 7 



Recently, we have made progress on this problem for the Vlasov equation; both in the 
approach to constructing the expansions and in the way the singular limit 7 — > + is treated 
and interpreted. |[23|| - [ 25[j An amplitude equation for an unstable mode in a one-dimensional 



collisionless plasma was derived for the dynamics on the two-dimensional unstable manifold 
of the equilibrium E . The essential difference from previous work lies in the choice of 
unperturbed state. Earlier theories assumed an equilibrium with a neutrally stable mode 
and obtained ill-defined expansion coefficients. |13[] This can be avoided by taking the weakly 



unstable equilibrium as the unperturbed state; a choice that naturally leads one to work 
with the unstable manifold. The mode eigenvalue A = 7 — iu can be complex (beam-plasma) 
or real (two-stream); in either case the equations for the amplitude A(t) = p(t) e^ te ^\ 

p = lp + a lP 3 + a 2 p 5 + 0(p 7 ) (1) 
= lu + a[p 2 + a' 2 p 4 + 0(p 6 ) (2) 

lead to a one-dimensional problem for p because the spatial homogeneity of the equilibrium 
decouples the phase dynamics. As 7 — > + , the expansion coefficients diverge, 
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"3> "3 7 4j-i 



(3) 



but these divergences can be removed to all orders in the expansion by rescaling the mode 
amplitude: p(t) = / ~f 2 r( / ~ft). In this way one obtains asymptotic equations for r(r) that are 
well behaved as 7 — ► + , and moreover through Poisson's equation this rescaling implies 
that the electric field exhibits the trapping scaling. These initial results were obtained for a 
plasma of mobile electrons with infinitely massive ions providing a fixed neutralizing back- 



ground; consequently they contain no information regarding unstable ion acoustic modes, 
for example. 

In order to study the effects of ion dynamics, we have generalized the analysis to treat 
an unstable electrostatic mode in a multi-species one-dimensional plasma. This changes the 
problem in a qualitative way: now as 7 — > + , the single particle dynamics and collective 
motions of the ions occur on a fast time scale relative to I/7. We calculate explicit ex- 
pressions for the leading nonlinear coefficients a\ and a[; in addition we have determined 
the dominant singularities in the amplitude expansion to all orders. The results show a 
qualitatively different singularity structure from the limiting model (H) with fixed ions, and 
provide new predictions for the scaling of nonlinearly saturated modes. 

The theory is described for the simplest example of a neutral plasma with two species 
(s = e, i) which we refer to as "electrons" and "ions" although the results apply equally to 
collisionless electron-positron plasmas. Let n s = N s /L denote the average species density in 
a one-dimensional plasma of length L, and eq s denotes the charge per particle of species s. 
In convenient dimensionless variables, the Vlasov-Poisson system becomes 

/oo 
dvF {s) (4) 
-00 

where kS s > = (q s m e /m s ). We assume periodic boundary conditions and adopt the normal- 
ization 

[ L/2 dx I" dvF^(x,v,t) = ^L; (5) 
J-L/2 J-00 n e 

note that F^ is negative for electrons and positive for ions. Given a spatially homoge- 
neous equilibrium Fq S \v), this system determines an evolution equation for f( s \x,v,t) = 
F&(x,v,t) -F^ s \v) 

ft/W =£/(•)+ JV(/W) (6) 

where Cf^ = -vd x f^ - K^Ed v F^ s) and Af(f {s) ) = -K^Ed v f^'\ 
An unstable mode exists if the dielectric function 
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ek W =1 -wL dv — ~ — 



(Im z > 0) 



(7) 



has a root £o — Vp + i'y/k in the upper half plane (7 > 0). The root determines an eigenvalue 
A = —ikzo for L with a two- component eigenvector 



* = e 



(8) 



/ 



We assume that there is a single such mode and that it corresponds to a simple root of €k(z), 
i.e. e' k (zo) 7^ 0. The eigenvalue A can be real or complex depending on the equilibrium; the 
ion-acoustic instability corresponds to a complex eigenvalue. 

The amplitude of the unstable mode is the coefficient of \I/ in the expansion of / 



f(x, v, t) = [A(t)ty(x, v) + cc] + S(x, v, t); 



(9) 



here / denotes the two-component field / = (f^ e \f^)- This decomposition allows the 
dynamics of the mode amplitude A(t) and the remaining modes S(x,v,t) to be separated 



A = XA+(^,Af(f)) 
d t S = CS + M(f) - ti% Af(f)) * + cc 



(10) 
(11) 



using an inner product defined for two-component fields B = (B^ e \ B^') and D = (D^ e \ D^) 
by (B,D) = J dx J dv [B^ e \x, v)*D^ e \x, v)+B®(x, v)*D i:i \x, v)], and the adjoint eigenfunc- 
tion \l/ for A*. In (|10])-(|ll]) the linear terms are decoupled, but nonlinear couplings between 
A and d t S remain. 

The amplitude equation for A follows when we express the time dependence of S in terms 
of A. 



S(x, v, t) = H(x, v, A{t), A*(t)). 



(12) 



As we have discussed elsewhere, this step can be visualized as a restriction of the initial 
condition to the two-dimensional unstable manifold of the equilibrium. ]24j2^| Consistency 
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between the time dependence of S in (|i~2| ) and the evolution of 5" described by (|10|)-(|TT 
requires that H(x,v, A(t), A*(t)) satisfy 



Ad A H + A*d A *H 



= CH + N{f u )- (*,Mf))*+ {^,M{f u )Y^* (13) 
/=/« L J 



where f u (x,v) = [A^(x,v) + cc] + H(x,v, A, A*). For solutions of this form, the dynamics 
of A(t) (|10D yields an autonomous equation for A 

A = \A+{$,N{f u )) (14) 



provided H can be determined from ([13|) . The homogeneity of the equilibrium i*o forces 
this amplitude equation to have a simple form: A = Ap(\A\ 2 ) where p(|A| 2 ) must still be 
determined. In polar variables, A = pe~ ld , the system (|i~4[ ) separates 

p = pRe[p(p 2 )} P 9 = -lm[p(p 2 )} (15) 

yielding a one-dimensional flow for p(t); the essential problem is to study p(p 2 )- 

Our conclusions regarding the evolution of the wave are based on an analysis of the 
amplitude expansion for p, 

oo 

p(p 2 ) = Y.pjp 2j > (16) 

j=0 

and similar expansions for H(x, v , A, A*). By substituting A = Ap(p 2 ) into ( |H]) we obtain 
one set of relations between the coefficients of p and H: 

oo 

AY,Pif?* = \A + (*,N{r))\ (17) 
i=o 

the defining equation (|13|) for H provides a second set of relations. The expansion coefficients 
are calculated by solving (^) and (|T7| ) recursively; full details of this calculation and our 



analysis of the resulting recursion relations will be published elsewhere. 

The first coefficient p$ is simply the linear eigenvalue A; information about the nonlinear 
evolution is provided by the higher order coefficients, in particular the first nonlinear term 
is given by 
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Pi = — 

T 



k«(1 - ^) 2 )Im(a(z )) /l k«(1-k«V(*o) , , 

7 7 H 77T77 a + U W 



18) 



Ake' k (z ) '[A 4k%(z ) 

where a'(z) is the derivative of a(z) = / dvd v F^ / (v — z). If the first term in (|T8|) is non-zero, 



i.e. 



— K (i > 2 )Im(o;(^o)) t^O, (19) 



then in the limit of weak instability 7 — > + , there is a 7 -4 singularity in p\. When ( |19"D 
holds, the trapping scaling introduced previously for p no longer removes the singularity 
in pi, and a new scaling pit) = 7 5//2 r(7t) is required to obtain a nonsingular leading order 
term. 

Beyond this leading order calculation, we have proved that the coefficients in ([16]) have 
the asymptotic behavior pj ~ to all orders j > 1. This general result implies that 

when the equations ([15]) are rewritten in the rescaled amplitude r the singularities will be 
cancelled to all orders in r. [2B| Thus given the generic condition ( |T9"D the theory predicts 



the asymptotic scaling for the mode electric field is \Ef.\ ~ 7 5 ^ 2 as 7 — > + . 

Our prediction must be re-examined if the genericity condition ( |T9"D fails and the leading 



term in p\ vanishes. The three factors in (TBI) indicate three exceptional circumstances when 



this can happen: «W = 0, Im(a(zo)) — or k^ 2 = 1. The first circumstance corresponds to 
the limit of fixed ions, m e /mj — > 0, and the cubic term (0) reduces to our previous result 
(0). With fixed ions, the terms of the amplitude expansion are known to be nonsingular to 
all orders once the amplitude has been scaled with 7 to balance the divergence in pi, i.e. 
p{t)=l 2 r{ 1 t). m 

The second circumstances arises naturally when the equilibria under consideration are 
reflection symmetric Fq S \—v) = Fq(v). With such a reflection symmetry, one can find 
pure imaginary roots Zq = —z (hence real eigenvalues) and then a(z ) is forced to be real; 
an example is provided by a reflection-symmetric two-stream instability. In this case, the 
trapping scaling for p cancels the singularity in p\ and moreover can be proved to yield an 



amplitude expansion which is nonsingular to all orders. |26| As in the fixed ion case, this 



result predicts the familiar trapping scaling for the mode electric field |i£fc(£)| ~ 7 2 . 

The third exceptional circumstance k^ 2 = 1 requires (gjm e /mj) 2 = 1 which corresponds 
to an electron-positron plasma (qi — 1, rrii — m e ). In this case the singularity structure 
of the expansion is more complicated. The cubic coefficient has a 7~ 3 singularity which 
suggests trapping scaling for p, but the divergence of fifth order coefficient turns out to 



be P2 ~ 7 -8 which is not removed by trapping scaling. |26j This fifth order singularity is 
cancelled if we set pit) = 7 9//4 r(7t). However, inspection of the higher order terms in the 
expansion shows the divergence structure of pj ~ 7~( 5j-3 ) and cancelling these singularities 
to all orders again requires the generic scaling pit) = 7 5//2 r(7t) even though to any finite 
order pj a smaller exponent would suffice. fl26f 

With the notable exception of this third example, our conclusions are easily summa- 
rized: the scaling required to obtain a nonsingular expansion is correctly predicted by the 
divergence found in the cubic coefficient p\. Generically, this singularity dictates a scaling 
by 7 5//2 , but this is replaced by 7 2 in the limit of infinite ion inertia or for instabilities in 
reflection-symmetric systems due to real eigenvalues. In the generic case we can estimate 
the range of growth rates where the new 7 5//2 scaling is visible by determining the range 
in 7 where the 7~ 4 divergence dominates p\. In general, this will depend on the specific 
parameters of the equilibrium, and we give two illustrative examples: an unstable plasma 
wave driven by an electron beam and an ion-acoustic instability. For simplicity, we have 
evaluated the dispersion relation and the cubic coefficient using Lorentzian distributions for 
the two species, 



l(e) _ -a e n p /ir a b n b /7T (i) _ K^oti/ir 

[v — u e ) 1 + aj iv — UbY + a^ v 2 + af 



with n p + rib = rii = !■ 

The plasma wave instability is calculated for an equal density cold beam with four mass 
ratios; the linear stability boundaries are shown in Fig. 1. The densities in this example 
are 2n p = 2n b = n» = 1, and u e = 0. The variation of 7 3 Repi as 7 — > + is shown in Fig. 
2. For each mass ratio in Fig. [I], we fix k and vary Ub] the chosen values are k = 0.75 for 
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m e /mi = 0.5 and k = 0.5 for m e /rrii = 0.1,0.01, and 0.001. As m e /rrii decreases, the fixed 
ion result 7 3 Rep 1 — > — 1/4 holds down to smaller and smaller growth rates. Interestingly, 
the effect of the mobile ions is to shift the asymptotic sign of Repi from negative to positive. 
Since the higher order terms in the amplitude expansion cannot be neglected, this does not 
automatically imply the onset of subcritical bifurcation but may nevertheless be significant. 

The ion-acoustic instability is shown in Figs. 3 also for four mass ratios. In this case 
rib = and there is only a drifting electron population with n p — rii — 1. In Fig. 4, 
for each mass ratio, we fix k and vary u e ; the chosen values are k = 1.5, 1.2,0.75, and 0.4 
corresponding to m e /mi = 0.5, 0.1, 0.01, and 0.001, respectively. Here the effect of decreasing 
m e /mi is to increase the range of the generic scaling as measured in units of the ion plasma 
frequency uf = Ane 2 rii/mi. 

The qualitative result of the scaling \E k \ ~ 7 5 / 2 is to reduce the electric field of the 
nonlinearly saturated wave compared to a wave characterized by trapping scaling. More 
work is needed to understand how the mobile ions manage to cut off the growth of the 
unstable mode at such small amplitudes. 
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FIG. 1. Beam-plasma stability boundary for a e = 1.0, oti = 1.0 and = 0.1 at four 
mass ratios m e /mi = 0.001,0.01,0.1 and 0.5. The boundaries for m e /mi = 0.001 and 0.01 
coincide. 
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FIG. 2. Asymptotic behavior of 7 3 Repi for the beam-plasma instability of Fig. [l]. 
The growth rate 7 is measured in units of the electron plasma frequency. The diver- 



gence as 7 — > + indicates the regime predicted to show the generic scaling 
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FIG. 3. Ion-acoustic stability boundary for a e = 1.0, a« = 0.01 at four mass ratios 
rrie/mi = 0.001,0.01,0.1 and 0.5. 
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FIG. 4. Asymptotic behavior of 7 3 Re(pi) for the ion-acoustic instability of Fig. ^ 
The growth rate 7 is measuredin units of the ion plasma frequency. The divergence 
as 7 — » + , indicates the regime predicted to show the generic scaling \Ek\ ~ 7 s / 2 . 
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